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Abstract 

We develop a rigorous quantum mechanical theory for collisions of polyatomic molecular radicals 
with S-state atoms in the presence of an external magnetic field. The theory is based on a fully 
uncoupled space-fixed basis set representation of the multichannel scattering wavefunction. Explicit 
expressions are presented for the matrix elements of the scattering Hamiltonian for spin- 1/2 and 
spin-1 polyatomic molecular radicals interacting with structureless targets. The theory is applied 
to calculate the cross sections and thermal rate constants for spin relaxation in low-temperature 
collisions of the prototypical organic molecule methylene [CH2(X 3 £?i)] with He atoms. To this 
end, two highly accurate three-dimensional potential energy surfaces (PESs) of the He-CH2(X 3 i?i) 
complex are developed using the state-of-the-art CCSD(T) method and large basis sets. Both PESs 
exhibit shallow minima and are weakly anisotropic. Our calculations show that spin relaxation in 
collisions of CH2, CHD, and CD2 molecules with He atoms occurs at a much slower rate than 
elastic scattering over a large range of temperatures (1 [iK - IK) and magnetic fields (0.01 - 
1 T), suggesting excellent prospects for cryogenic helium buffer-gas cooling of ground-state ortho- 
CH2(X 3 i?i) molecules in a magnetic trap. Furthermore, we find that ortho-CH2 undergoes collision- 
induced spin relaxation much more slowly than para-CH2, which indicates that magnetic trapping 
can be used to separate nuclear spin isomers of open-shell polyatomic molecules. 
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I. INTRODUCTION 



Cooling and trapping polyatomic molecules is a promising research direction within the 
rapidly expanding field of cold molecular gases 1 . Because of their rich internal structure 
and coupled vibrational degrees of freedom, polyatomic molecules offer new research oppor- 
tunities at the interface between chemistry, physics, combustion 2 , atmospheric chemistry 3 , 
and astrochemistry 4 . Studies of intramolecular energy redistribution in complex organic 
molecules yield insight into a broad range of phenomena, ranging from biochemical reac- 
tions involved in visual transduction 5 to quantum decoherence in mesoscopic systems 6 . Un- 
like diatomic molecules, polyatomic molecules composed of four (or more) atoms can exist 
in the form of left and right-handed isomers and high-precision spectroscopy of cold molec- 
ular ensembles can provide insight into parity violation 7 and Hund's paradox 8 . Internal 
cooling increases the population of low-lying rovibrational states and hence the molecular 
response to external field perturbations, while cooling of external (translational) degrees of 
freedom enhances the molecule-field interaction time. High-precision spectroscopy of cold 
polyatomic molecules may thus enable ultrasensitive molecular detection and separation of 
different conformers 9 as well as the development of novel time and frequency standards 10 . 

While indirect cooling techniques such as photo-association and magneto-association are 
highly efficient at producing ultracold alkali-dimers 11 , they cannot be extended in a straight- 
forward manner to cool polyatomic molecules due to the lack of efficient laser cooling meth- 
ods for all but a handful of atomic species 1 . Direct cooling techniques such as sympathetic 
cooling 12-15 , velocity filtering 16 , rotating nozzle slowing 17 , and Stark deceleration 18 are free 
from these difficulties, as demonstrated by recent experiments on the production of cold en- 
sembles of fully deuterated ammonia 19 , benzonitrile 20 , formaldehyde 21 , and naphthalene 22 . 
Slow beams of polyatomic molecules thus produced have already been used, in combination 
with electromagnetic traps, to explore the dynamics of OH + ND 3 23 and Li + SF 6 24 colli- 
sions at low temperatures. However, quantum collision dynamics of polyatomic molecules 
in the presence of external electromagnetic fields remains poorly understood, necessitating 
the development of theoretical tools for ab initio modeling of these pioneering experiments 
as well as for the realization of efficient strategies for controlled molecular cooling 1 . 

Theoretical studies of polyatomic molecule collisions in the absence of external fields 
were pioneered by Green 25 , who developed a formal quantum theory of collisions involving 
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asymmetric tops, and carried out scattering calculations on various problems of astrophysical 
interest, including He + HCN, He + H2O, and He + NH3 25-27 . Many authors employed 
Green's formalism to carry out accurate scattering calculations of rotational cooling rates 
for interstellar molecules colliding with atomic and molecular species based on accurate ab 
initio interaction potentials 28 . Several theoretical studies have used a combination of high- 
level ab initio and quantum scattering calculations to explore the dynamics of collision- 
induced rovibrational relaxation of methylene (CH^a, X)) and methyl (CH3) radicals in 
room-temperature He buffer gas 29,31 , for which accurate rate measurements have recently 
become available 32 . 

A number of related theoretical studies have focused on field-free collisions of H2O, CO2, 
and benzene molecules with He atoms at temperatures below 1 K 33,34 . These calculations 
were motivated by the ongoing experiments on sympathetic cooling of optically decelerated 35 
polyatomic molecules with rare gas atoms 34 . In particular, evidence was reported for slow 
rotational energy transfer in cold collisions of highly rotationally excited CO2 molecules with 
He atoms 36 . A recent classical trajectory study of He-naphthalene scattering 37 has shown 
that the unexpected lack of cluster formation observed in buffer-gas cooling experiments on 
naphthalene 22 can be attributed to the short lifetime of the He-naphthalene complex. 

We have recently presented a rigorous quantum scattering methodology for numerical 
calculations of collisional properties of polyatomic molecular radicals with 5*-state atoms in 
the presence of an external magnetic field 38 . The method is based on the fully uncoupled 
space-fixed representation of the scattering wave function originally introduced by Volpi 
and Bohn 39 and Krems and Dalgarno 40 to investigate collisions of diatomic molecules in a 
magnetic field. In Ref. 38 we studied the quantum dynamics of spin relaxation in collisions of 
spin-1/2 and spin-1 polyatomic molecular radicals induced by collisions with He atoms. We 
found that many of these polyatomic molecules have favorable collisional properties (that 
is, large ratios of elastic to inelastic collisions rates) with He atoms, making them amenable 
to cryogenic buffer-gas cooling and possibly sympathetic cooling with S'-state atoms in a 
magnetic trap 38 . 

Here, we present a full account of the theory and computational methodology outlined in 
our previous work 38 . We begin with the derivation of coupled- channel quantum scattering 
equations for a polyatomic molecule colliding with an 5-state atom in the presence of an 
external magnetic field (Sec. IIB). In Sec. Ill, we use the newly developed scattering 
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methodology to address the question of collisional stability of methylene radicals in cold 
3 He gas in the presence of an external magnetic field. After a brief outline of the ab initio 
techniques employed to obtain the adiabatic PES for the He-CH 2 complex (Sec. IIIA), we 
present and discuss the main results for elastic and inelastic cross sections and thermal 
rate constants for spin relaxation in He + CH 2 collisions (Sec. IIIB). We then analyze the 
scaling properties of the inelastic cross sections and discuss the mechanism of spin relaxation 
of different nuclear spin isomers (Sec. IIIC) and isotopomers (Sec. HID) of CH 2 . Section 
IV summarizes the main results of this work. 

Atomic units are used throughout unless stated otherwise. 



II. THEORY 



In this section, we will outline our approach to the solution of the atom - polyatomic 
molecule collision problem in the presence of an external magnetic field. We first describe 
our ab initio calculations of the potential energy surfaces (PESs) of the He-CH 2 complex. In 
Sections IIB and IIC, we generalize this work to collisions of polyatomic molecules (asymmet- 
ric tops) with S'-state atoms. We present the derivation of the matrix elements required to 
parametrize the close-coupling (CC) scattering equations for collisions of polyatomic molec- 
ular radicals with S*-state atoms. Sec. II C presents a brief description of our numerical 
calculations. 



A. Hamiltonian 

The Hamiltonian for an open-shell polyatomic molecule colliding with a 1 S'o atom in the 
presence of an external magnetic field may be written as 

where /i is the reduced mass of the collision complex, R = RR is the atom-molecule separa- 
tion vector of length R, I 2 is the orbital angular momentum for the collision, V(R, Cl) is the 
anisotropic potential describing the atom-molecule interaction, and H mo \ is the asymptotic 
Hamiltonian describing the internal structure of the molecule and its interaction with exter- 
nal fields. In this work, we focus on collisions at low incident kinetic energies (10 -4 - 100 K), 
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much lower than the typical energy scale for intramolecular electronic and vibrational exci- 
tations (typically, hundreds and thousands of Kelvin) so we can assume that the molecule 
remains in its ground electronic and vibrational states during the collision. Thus we can 
use the following effective Hamiltonian 42-44 to describe the rotational and fine structure of 
a polyatomic molecular radical like CH 2 45 : 

-f^mol = -f^rot + H C( ± + H ST + H ss + H ext . (2) 

The rotational part of this Hamiltonian is that of a rigid asymmetric top with rotational 
constants A, B, and C: H rot = AN 2 + CN£ + BN 2 , where the operators N a (a = a,b,c) 
yield the molecule-fixed (MF) components of rotational angular momentum N. The fine- 
structure Hamiltonians H ST and H ss describe the spin-rotation and spin-spin interactions in 
polyatomic molecules with non-zero total spin S (radicals). The term H c & describes the 
centrifugal distortion of rotational levels and the term H ext describes the interaction of the 
molecular spin with an external magnetic field. 

Here, we choose the MF axes to correspond to the symmetry axis of the molecule, as 
illustrated in Fig. 1. This choice is the same as adopted by Hutson in his theoretical study 
of Ar + H2O collisions but differs from the recent work of Dagdigian and co-workers, where 
the principal axis of the triatomic is used to define the z axis 29 ' 30 . The CH 2 molecule lies in 
the xz plane, with the z-axis bisecting the HCH angle, and the y-axis perpendicular to the 
molecular plane. With this choice, the He-CH2 interaction potential is symmetric under the 
reflection in the xz plane. The rotational part of the effective Hamiltonian (2) becomes 

Hrot = AN 2 + CN 2 + BN 2 , (3) 

and the centrifugal distortion Hamiltonian is (up to fourth order) 42 

H cd = -A N N* - A NK N 2 N 2 X - A X A> 4 - 2«5 Jv iV 2 (iV 2 - N 2 y ). (4) 

The spin-rotation interaction is given by 42 ' 48 

^sr = lxN x S x + lyNySy + i z N z S z , (5) 

where S is the electron spin, 7 a are the spin-rotation constants. Eq. (5) gives the spin- 
rotation interaction in terms of the MF angular momentum operators N and S. The 
corresponding expression in terms of SF angular momentum operators is derived in Appendix 
A. 
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The spin-spin interaction may be written as 

2 r 

# SS =E \{E + D)[Vl 2 (h) + Vl_ 2 {h)\ + ?*^z£ (n) [SI ® (6) 

where 2^ )P ,(fi) are the Wigner .D-functions, f2 are the Euler angles which specify the ori- 
entation of MF axes in the space-fixed (SF) coordinate frame, and D and E are the zero- 
field-splitting parameters 42 ' 48 . Equation (6) may be derived from the standard textbook 
expression for the spin-spin interaction in the MF frame 42 ' 48 ' 49 as described in Appendix A. 
Note that the spin-spin interaction is only different from zero in molecular radicals bearing 
more than one unpaired electron (S — 1). 

The interaction of the polyatomic molecule with an external magnetic field is described 
by the Hamiltonian 42 

#ext = /iO-B • g s ■ S, (7) 

where B is the magnetic field vector, /i is the Bohr magneton, and g s is the g-tensor. 
Following Ref. 42 , we take g s = 21 where I is the unit matrix, so Eq. (7) reduces to a form 
commonly used for diatomic 3 S molecules 

Hext = 2/io-BSz, (8) 

where Sz yields the projection of S on the magnetic field axis. The rotational and nuclear 
spin Zeeman interactions are much smaller than the H ex t term, and we neglect them in this 
work. Table II lists the molecular constants used in this work to parametrize the Hamiltonian 
(2) for CH 2 (X). 



B. Coupled-channel equations, S-matrix elements, and cross sections 

To solve the quantum scattering problem for the Hamiltonian (1), we expand the wave 
function of the collision complex in direct products of internal basis functions \a) and partial 
waves of relative motion \lrrii) in the space-fixed (SF) coordinate frame 

l*> = 5 E F almi (R)\a)\l mi ), (9) 

a,l,mi 

where the basis functions |a) = \NM^Kn)\SMs) describe the rotational and spin de- 
grees of freedom of the polyatomic molecule, \NM^Kn) = a/ (2N +1)/ 87r 2 "Z>^ K (Cl) is a 
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symmetric-top eigenfunction, Mn (-Kjv) is the projection of N on the SF (MF) quantization 
axes, and Ms is the SF projection of S. We note that the projection M = M^- + Ms + mi of 
the total angular momentum on the SF quantization axis is rigorously conserved for collisions 
in a magnetic field. As in the atom-diatom case considered previously 40 this conservation 
can be exploited to factorize the collision problem into smaller problems which can be solved 
independently for each M. 

Using the expansion (9) in combination with the time-independent Schrodinger equation 
for Hamiltonian (1) leads to a system of close-coupling (CC) equations 

F almi (R) = 2fi («|(/m i |y(/?,^0) + if mo i|a , )l^)^' m ;( J R)- 

a 1 ,V ,m'i 

(10) 

The matrix elements of the asymptotic Hamiltonian H mo \ and of the atom-molecule inter- 
action potential V(R, 9, 0) are evaluated as described below. As in the atom-diatom case 40 , 
the matrix of the asymptotic Hamiltonian is not diagonal in the fully uncoupled SF repre- 
sentation. In order to properly apply the scattering boundary conditions, we introduce an 
asymptotic basis 40 

\j)\lmj) = \lmi) y^C a7 |a), (11) 

a 

which has the property 

(Tl-^mollY) = <V/e 7 , (12) 
where e 7 are the asymptotic eigenvalues, which define scattering channels in the presence 
of an external magnetic field. The coefficients C ai in Eq. (11) are independent of I and 
mi and form the matrix C, which satisfies C T H mo iC = E, where E is the diagonal matrix 
of eigenvalues e 7 . The asymptotic behavior of the solutions to CC equations (10) in the 
asymptotic basis 40 defines the elements of the scattering S"-matrix 

FM m ,{R) ~ 6^6 u ,6 mm/ e-^ R -^ - ^ S^ m ^ R ~^ {R -+ oo), (13) 

as in the atom-diatom case 40 . The cross sections for the collision-induced transition 7—7-7' 
in the polyatomic molecule may be expressed through the S'-matrix elements as 

CT^y(E c ) = T2 ^2 ^2 l^rr'^f'^mj - ^mi;Yi'm[( E c)\ 2 , (14) 

7 M l,m l V,m' l 

where k 2 ^ = 2fi(E — e 7 ) is the wavevector for channel 7, E = Ec + e 7 is the total energy, and 
Eq is the collision energy. 
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d 1 

dR 2 



1(1 + 1) 

R 2 



2fiE 



C. Matrix elements in the fully uncoupled space-fixed basis 



The matrix elements entering the CC equations (10) can be expressed analytically via 
molecular constants and expansion coefficients of the atom-molecule interaction potential. 
In this section, we present the expressions for these matrix elements in the fully uncoupled 
SF basis of Eq. (9). 

We start with the asymptotic Hamiltonian of the polyatomic molecule (2). First, we note 
that since the asymptotic Hamiltonian (2) is independent of the atom-molecule orientation, 
its matrix elements are independent of I and m^: 

(a^lm^H^a'^l'm'i) = 5^5 m;m j (a|# mo i|a'). (15) 

The matrix element on the right-hand side can be written as the sum of rotational (3), 
centrifugal distortion (4), fine-structure (5), (6) and magnetic field-dependent terms (8) 
given by Eq. (2). 

In order to evaluate the matrix elements of H TOt , we express the MF angular momentum 
operators N x and N y in terms of the ladder operators N± and use the well-known expressions 
for their matrix elements 46 to obtain 

(NM N K N \(SM S \ H TOt | N'MffK'ff) | S M' s ) 

'1 



&M S M' S &NN' &M N M' N 



-{A- C)X ± (N, K' N )X ± (N, K' N ± 1)5 Kn , k , n±2 



+ \{A + C)[N(N + 1) - K 2 n ]5 KnK > n + BK 2 n 5 KnK , 



(16) 



where 

X ± (N, K N ) = [N(N + 1) - K N (K N ± l)] 1 / 2 . (17) 

The matrix elements of the centrifugal distortion term (4) can be evaluated following the 
same procedure. Explicit expressions for the matrix elements of H c ^ are presented in Ap- 
pendix B. 

A triatomic molecule with two identical H nuclei (like CH 2 ) can exist in two different 
nuclear spin modifications, ortho (o) and para (p). Nuclear spin selection rules strongly 
suppress collision-induced transitions between the energy levels of different nuclear spin 
isomers 48 ' 50 . Because the Hamiltonian is invariant to interchange of the two H atoms, o- 
CH 2 and p-CH 2 can be considered separate molecular species 26 . We obtain the energy 



level spectrum of each nuclear spin isomer by diagonalizing the Hamiltonian (2) using a 
restricted basis of symmetric-top eigenfunctions \N M^Kn) with (—1)^ = 1 (for 0-CH2) 
and (-1)*" = -1 (forp-CH 2 ). 

In order to evaluate the matrix elements of the spin-rotation interaction (5), we rewrite 
it in the form 



Hta 



S+i 

S-i + j z N S Q , 



(18) 



where N q and S q are the spherical tensor components of N and S in the MF frame. Rewriting 
the MF spin operators in terms of their SF counterparts 46 (see Eq. (35) in Appendix A), 
evaluating the integral over three Wigner Z)-functions, and rearranging the result, we obtain 



(NM N K N \(SM s \H sr \N'M' N K' N )\SM' s ) = [{2N + 1)(2N' + 1)] 1/2 [(2S + 1)S(S + 1)] 1/2 

S 1 S \ f N 1 N' 
-M s q M'sj \M N q -M' r 



x(-l) 



M'-K'+S- 



M *M 



NK N ;N'K' N ^ 



(19) 



iV 



where 



■MnK n ;N'K' 



2V2 . 



2V2 



\_{N',K' N ){ lx - ly ) 
+\ + (N\K' N )( lx + ly ) 
\+{N',K' N ){ lx - ly ) 
+\4N>,K' N )( lx + % ) 



N 1 N' 
K N 1 -K' N + 1 

N 1 N' 
K N 1 -K' N - 1 



(20) 




N 



N 1 N' 
K N -K' 



N 



and the quantities in parentheses denote 3-j symbols. An alternative expression for the 
matrix elements of the spin-rotation interaction may be obtained by directly evaluating the 
matrix elements of the tensor product in Eq. (5) in the SF frame (see Appendix A). The 
resulting expression is rather cumbersome, and we do not present it here. We used this 
expression to cross-check Eq. (20) and verified that both expressions yield the same results. 
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The matrix elements of the spin-spin interaction (6) in the fully uncoupled basis factorize 
into products of matrix elements of the form 

{NM N K N \Vl 2 (Q)\N'M' N K' N ) (SM S \[S ® S}® | SM' S ) . (21) 

The first matrix element reduces to a product of two 3-j symbols 46 and the second matrix 
element is given by Eq. (19) of Ref. 40 . Combining these results, we obtain (neglecting 
constant ^-dependent terms) 

(NM N K N \{SM S \H SS \N'M' N K' N }\SM' S } = [(2N + 1)(2N' + i)^/^-1)M' n ~k' n +s-m s 



\<E + D) 




N 2 N' 
K N -2 -K' N 



3E-D 




V(R, 6,4>) = W V Xll (R)——[C x ,(9, 0) + {-iyC x ^{9, 0)], (23) 



1 1 2 JV 2 JV 

x V5[(2S + 1)S(S + 1)] I \y\ 

[ s s s J Y V Mn q - M * 

Substituting = in Eq. (22), we recover the expression for the matrix element of the 
spin-spin interaction derived previously for linear 3 S molecules 40 . 

In order to evaluate the matrix elements of the atom-molecule interaction potential 
in the fully uncoupled SF basis (9), we expand the potential in renormalized spherical 
harmonics 25,41 

EE- 

A=0 m>0 

where 9 and are the spherical polar angles of R in the MF frame (see Fig. 1) and 
C*Aft(^)0) = [47r/(2A + \)] l l 2 Y\^{6, 0). The expansion (23) is valid for any polyatomic 
molecule interacting with an S'-state atom provided that V(R,9,(f)) is symmetric under 
reflection in the xz plane, i.e., V(R,6,4>) = V(R, 8, — 0). In our case, because of the C-i v 
symmetry of CH2, the potential is also symmetric under reflection in the yz plane, so we 
have V(R, 6, 0) = V(R, 6, it — 0). This condition implies that the expansion coefficients Vx^ 
vanish for odd \x. We note that this property is a consequence of our using a specific MF 
system shown in Fig. 1. As mentioned above, Dagdigian and co-workers chose their MF 
z-axis to be perpendicular to the C2 axis of CH2 29 - 31 . As a result, their coefficients are 
not the same as defined in this work, and they do not vanish for odd /i 29 ~ 31 . 

Because our SF basis functions do not depend explicitly on 9 and 0, we transform the 
expansion (23) to the SF frame. Transforming the spherical harmonics to the SF system, 
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we obtain 25 

A=0 v>0 M ° m x 

While this expansion is convenient for obtaining analytical expressions for matrix elements, 
it involves five angular variables, and is thus not well suited for computational purposes. A 
convenient way of evaluating the coefficients is outlined in Appendix C. 

Making use of standard expressions for the integrals involving products of three D- 
functions and spherical harmonics 46 , we obtain the final result for the matrix element of 
the interaction potential 



(NM N K N \(SM s \(l mi \V(R, Q)\N'M' N K' N )\SM' S ) |/'m|) = 

5 MsiM ,[(2iV + l)(2iV / + l)(2/ + l)(2/ , + l)] 1 / 2 x (-1)^-^+^ 

N X N' 
M N m x —M' N 

N X N' 
K N /i -K' N 



x 

X 



A=0 n>0 ^ m A \ 




(25) 



Setting Kn = reduces this expression to the familiar Legendre expansion of the atom- 
diatomic molecule interaction potential 40,51 . As in the atom-diatom case 40 , the interaction 
potential V(R, Cl) does not couple basis functions with different Ms- 

The coefficients Vx^R) were evaluated as described in Appendix C using a 15 x 15 direct- 
product Gauss-Legendre quadrature grid in 9 and 0. The expansion (25) was truncated at 
Amax = 7 and we verified that increasing A max does not affect the numerical results. The 
system of CC equations (10) was solved numerically using the log-derivative algorithm 52 on 
a grid of R extending from -R m j n = 3.0ao to i? max = 40.04a with a grid step of 0.04ao for 
collision energies larger than 0.1 cm -1 . For ultralow collision energies (Ec < 0.1 cm -1 ), 
-Rmax was increased to 80.04a (for 10~ 4 cur 1 < E c < 0.1 cm" 1 ) and 160.04a (10" 6 
cm -1 < E c < 10" 4 cm" 1 ). 

To obtain the cross sections converged to within <10%, the coupled-channel basis set 
included 5 rotational states of CH2 augmented with 5 (for Ec < 1 cm" 1 ) or 6 (for Eq > 1 
cm" 1 ) partial waves. The same basis set parameters were used for all CH2 isotopomers. 
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III. RESULTS 



A. Ab initio PESs 

In this work, we use two different ab initio potential energy surfaces (PESs) for the He- 
CH2 van der Waals complex: PES A and PES B. The former (PES A) was utilized in our 
previous theoretical work on He + CH2 scattering in a magnetic field 38 . A more recent, 
high-quality PES B developed by Ma et al. 31 has been used in quantum calculations of 
rotational relaxation in He + CH 2 collisions at room temperature. 

PES A was calculated using a UCCSD(T) method (coupled-cluster including single and 
double excitations along with a perturbative correction for triple excitations) 54 and a large 
cc-pVTZ basis set 55 . For fitting purposes, the PES is described by an analytic function that 
sums over six pairs of interaction potentials, 53 

^He-CH 2 = Ys{ A « e ~ baRa ~ fdb«R«)^\, (26) 

with the damping function 

n l, 

/„(*) = 1 - X w (27) 

fc=0 

The six pairs refer to a = {He — C, He — H a , He — H^, He — Gi, He — G2, He — G3}, where 
Gj are three pseudo-atoms. Their MF coordinates relative to C are fixed at (0,1,0) a for Gi, 
(0,-1,0) a for G 2 , and (0,0,1) a for G 3 . The pseudo-atoms are used to take the many-body 
interaction potentials into account. 

The parameters in Eq. (26) are determined by fitting to 985 CCSD(T)/cc-pVTZ points 
spanning the range of R = [3.0, 7.5] a , 9 = [3°, 179°], and <p = [5°, 175°] with a fixed 
geometry of CH 2 at its triplet equilibrium structure. The potential energy surface is well 
fitted to a root-mean-square error of 0.8 cm -1 . The parameters obtained are given in Table I. 

Figures 2-4 (left panels) show a typical long-range potential energy surface for the 
He— CH 2 complex, where the CH 2 molecule is fixed at its CCSD(T) optimized geometries, 
Ren = 2.0404 ao and #hch = 133.11°. It displays two equivalent global minima located 
at i?Hc-CH 2 — 7.68 ao with a well depth of 9.7cm _1 . Owing to the shallow minima, the 
interaction potentials between the two collision partners are weakly anisotropic. This feature 
of the potential energy surface is also consistent with the small permanent electric dipole 
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moment of 0.2293 atomic units (0.5828 D) calculated with CCSD(T)/cc-pVTZ and nearly 
isotropic quadrupole moments [(6.0344, 4.3788, 5.6653) atomic units] of CH2. 

A detailed description of the ah initio approach and the topology of PES B has been 
presented elsewhere 31 . Spin-restricted coupled-cluster [RCCSD(T)] calculations 56,57 were 
carried out on a 4-dimensional grid consisting of the position (in spherical polar coordinates) 
of the He atom with respect to the center-or-mass of the CH 2 (AT) molecule for a grid of values 
of the CH2 bond angle. An aug-cc-pvqz 58 basis set was used with a set of bond functions 59 ' 60 
added at the mid-point of R, the vector pointing to the He atom from the center of mass 
of CH 2 . The 3-dimensional PES used in the scattering calculations was determined by 
averaging the 4-dimensional potential over the bending angle weighted by the square of the 
v b = bending vibrational wave function. 

Figures 2-4 show that PES A is considerably less deep, more repulsive at short range, 
and less anisotropic than PES B. The increased depth of the potential arises from a better 
recovery of the correlation energy, which is known to be responsible for dispersion interac- 
tions, due to the inclusion of diffuse atom-centered functions as well as bond functions 62 for 
a weakly-bound complex. The calculations on which PES A is based are deficient in that the 
van der Waals attraction is significantly less well described. We also note that because the 
van der Waals well is deeper in PES B, the repulsive core of the PES - the angle-dependent 
hard sphere radius - is significantly smaller in extent than the predictions of PES A. As will 
be shown below (Sec. IIIB), the significant difference in the two PESs has a strong effect on 
both the elastic and inelastic He + CH 2 (X) collisions at low temperatures. 

B. Low-temperature He + 0-CH2 collisions in a magnetic field 

Figure 5 shows the Zeeman energy levels of ortho and para nuclear spin isomers of CH 2 . 
The ground rotational state of o-CH 2 (Ooo) is split by magnetic fields into three Zeeman 
sub levels with Ms = —1, 0, and 1. Here, we adopt the conventional Nk k p notation for the 
rotational states of an asymmetric top, where K (K p ) are the projections of N on the MF 
quantization axis in the oblate (or prolate) symmetric top limits 67 . The ground rotational 
state of p-CH 2 (loi) nas N = l 42 - 47 and splits into nine Zeeman sublevels correlating with 
the j = 0, j = 1, and j = 2 fine-structure levels in the limit of zero magnetic field, where 
j = N + S is the total angular momentum of CH 2 excluding nuclear spin. Here, we are 
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interested in collisions of CH2 molecules in their magnetically trappable, maximally spin- 
stretched Zeeman levels Ms = +1 (for 0-CH2) and rrij = +1 (for p-CH 2 ), where mj = 
Ms + Mtv is the projection of j on the magnetic field axis. It is clear from Fig. 5 that the 
Zeeman energy level patterns of o-CH 2 and p-CH 2 are very different. As will be shown below, 
this difference leads to dramatic variations in the low-temperature collisional properties of 
the spin isomers. In this section, we focus on collisions of 0-CH2 molecules and omit the 
prefix o for brevity. 

Figure 6 shows the elastic and inelastic cross sections for 3 He + CH 2 (0 o, M s = +1) 
collisions 65 as functions of collision energy calculated for B = 0.1 T, the field strength typi- 
cally used in buffer-gas cooling and magnetic trapping experiments with 3 S molecules 14,15 . 
We calculate the inelastic spin relaxation cross section by summing the state-to-state cross 
sections for a given incident collision channel (dashed lines in Fig. 5) over all energetically 
accessible Zeeman levels of CH 2 (full lines in Fig. 5). Both the elastic and inelastic cross 
sections are very sensitive to the interaction PESs: the s-wave limit of the elastic cross sec- 
tion calculated using PES A is ~10 times larger than the limit for PES B. Overall, because 
calculations on which PES A is based predict a larger hard-sphere radius, we anticipate that 
the s-wave elastic cross section, which is a measure of the size of the system, will be larger. 

The cross sections calculated using PES A are dominated by a pronounced resonance 
feature at Eq ~ 0.1 K with a very broad shoulder extending down to Eq ~ 10 -4 K, whereas 
the results for PES B display two resonance peaks above Eq = 0.1 K. The presence of the 
broad resonance leads to a large difference between the thermally averaged ratios of elastic 
to inelastic collision rates calculated with potentials A and B (see the inset of Fig. 6). Since 
PES B has a more pronounced long-range attractive component, we anticipate that any 
elastic resonances at low energy predicted by the two PES's will be significantly different. 

Figure 7 shows the cross sections for inelastic spin relaxation in He + CH 2 collisions as 
functions of collision energy calculated for several values of the magnetic field. At very low 
collision energies, we observe the s-wave threshold behavior o"i ne i ~ l/\fEc typical of two- 
body inelastic processes in the limit of zero collision energy 64 . The elastic cross section is 
independent of the field and becomes a constant at Eq < 100 mK. At moderately low collision 
energies (above 0.5 K), the cross sections exhibit resonant oscillations due to the excited 
rotational states (/ > 0) of the collision complex (shape resonances) which decay by tunneling 
through the centrifugal barrier 15 ' 70 . A crossover between the single s-wave and multiple 
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partial-wave regimes occurs at a collision energy of 1 mK at low magnetic fields (B = 0.01 
T) and shifts to higher collision energies with increasing field. A similar trend manifests itself 
in cold collisions of 3 E molecules 40 . Figure 7 shows that the magnitude of the inelastic cross 
section in the s-wave limit increases dramatically with the field strength. As in previous 
theoretical work on collisions of E-state molecules, this can be explained by the following 
argument 39,40 . The height of the centrifugal barrier in the outgoing (Ms = 0, —1) collision 
channels becomes smaller with increasing field, making it easier for collision products to 
escape over the barrier, thereby enhancing the inelastic collision rate 39,40 ' 63 . 

In the inset of Fig. 7, we plot the key figure of merit for sympathetic cooling experiments 
- the ratio of the thermal rate constants for elastic scattering and inelastic scattering (spin 
relaxation) in He-CH 2 collisions: 

■< (T) = £S>- (28) 

The rates in Eq. (28) are obtained as functions of temperature by averaging the cross sections 
shown in Fig. 7 over a Maxwell-Boltzmann velocity distribution. We observe that 7 > 1000 
over the range of magnetic fields typically used in cryogenic cooling and magnetic trapping 
experiments (0.01 - IT). An empirical rule states that 7 > 100 is required to ensure efficient 
cryogenic He buffer-gas cooling of a molecule 12 , so the results shown in Fig. 7 indicate that 
p-CH 2 is an ideal candidate for cryogenic cooling and magnetic trapping experiments 12 . As 
shown in the inset of Fig. 6, the temperature profiles of 7 calculated for PESs A and B differ 
substantially, with 7a ^> 7b at T < 0.1 K due to the broad scattering resonance that occurs 
on PES A (see the inset of Fig. 6). We emphasize, however, that despite these quantitative 
differences, both PES A and PES B calculations predict j(T) > 100, so our main qualitative 
conclusion as to the suitability of CH 2 for sympathetic cooling experiments is independent 
of the choice of the PES. 



C. Spin relaxation of nuclear spin isomers 

So far we have been focusing on collisions involving a single nuclear spin isomer of CH 2 
(ortho-CH.2) . However, as pointed out in Sec. II, the rotational and Zeeman structures of o- 
and p-CH 2 are very different, so it is of interest to compare the collisional properties of the 
nuclear spin isomers. Figure 8 shows the elastic and inelastic cross sections for o-CH 2 and 
p-CH 2 in He as functions of collision energy. The inelastic cross section calculated for p-CH 2 
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is about four orders of magnitude larger than that for 0-CH2 over the range of collision 
energies from 10~ 5 to 0.1 K. In the multiple partial wave regime (Ec > 0.1 K) the difference 
becomes somewhat smaller, but never falls below three orders of magnitude. 

In order to understand this remarkable disparity in low-temperature collisional properties 
of different nuclear spin isomers of CH 2 , we note that in our calculations, the incident 
collision channels for both ortho and para-CH.2 are the maximally spin-stretched Zeeman 
states (shown in Fig. 5 by dashed lines). These are the Zeeman states in which the molecules 
reside when confined in a permanent magnetic trap 1,12 ' 63 . However, the rotational manifolds 
to which these Zeeman states belong are different for ortho and para nuclear spin isomers: 
As shown in the left panel of Fig. 5, the lowest maximally spin-stretched Zeeman state 
of 0-CH 2 belongs to the Ooo rotational manifold, while that for p-CH 2 belongs to the l i 
rotational manifold 47 . 

The spin-spin interaction (6) couples the ground Ooo an d the second excited 2 02 rotational 
manifolds of o-CH 2 . As a result, both N and Ms lose the status of good quantum numbers, 
and the wavefunction of the lowest magnetically trappable Zeeman state becomes 

|0 00 , Ms = 1) = a |0 00 , M s = 1) + O \2 O2 , M s = 0) + To|2 02 , Ms = -1), (29) 

where the tildes over N and Ms indicate the approximate nature of these quantum numbers. 
A first-order perturbation theory estimate gives 

o = (Ooo,M^ = 0|(ll|ff ss |2 02 ,M;,)|10) ^ _D_ 

AE02 ^Eq2 

where the basis functions are the products of asymmetric-top eigenfunctions \Nk k p Mn) 
times the spin functions \SMs), and we have neglected the dependence of the matrix element 
of the spin-spin interaction (6) on the zero-field-splitting parameter E assuming E <C D, 
which is justified for CH 2 (A > ) where E/D = 0.051. In Eq. (30) 67 

AE 02 = 2(A + B + C) - 2[{B - C) 2 + {A - C){A - B)] 1/2 (31) 

is the energy gap between the ground and the second excited rotational states. In 0-CH2 the 
gap (47 cm -1 ) is much larger than the spin- spin coupling D = 0.78 cm" 1 , so the admixture 
of the 2 2 rotational state into the ground state is on the order of a few percent (depending 
on M' N ). 

In the first Born approximation, the cross section for spin relaxation out of the lowest 
magnetically trappable state of CH 2 (|0 o,Ms = 1)) is proportional to the square of the 
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matrix element between the molecular spin states involved in the transition. In particular, 
the cross section for relaxation from the Ms = +1 state in the lowest rotational manifold 
of o-CH 2 may be written as 39 ctms^m' s ~ Iwooj M s \V(R 1 6>, 0)|O O o, M' s )\ 2 . Using the basis 
functions given by Eq. (29), we find for the Ms = 1 — >■ M' s = — 1 transition to the leading 
order 

<7i-o~0o x |(2 o2 |n^^0)|Ooo)| 2 ~ (^-) \{2o2\V(R,9,<f>)\0 00 )\ 2 . (32) 

This result shows that spin relaxation in spin-1 molecular radicals is a first-order process 
mediated by the combined action of the spin-spin interaction and the anisotropy of the 
interaction potential (the second term in Eq. 32). Note that both these interactions are 
equally important: without the spin-spin interaction there would be no admixture of the 
spin-down state in the initial spin- up state (D = implies /3q = in Eq. 30), and without 
the potential anisotropy there would be no coupling between the different N states. 

We further note that the expression (32) contains the factor D 2 /AE 2 2 <C 1 for light 
molecules with large rotational splittings (31). Thus, collision- induced spin relaxation in such 
molecules should be slower than rotational relaxation, which is induced by the anisotropic 
part of the interaction potential without the small factor D 2 /AE 2 2 . These considerations 
explain the small magnitude of the calculated inelastic cross sections for o-CH 2 shown in 
Fig. 7. A closely related result for collision-induced spin relaxation in linear 3 S molecules 
was first derived theoretically by Krems and co-workers 40 and then observed experimentally 
using He buffer-gas cooled NH( 3 S) molecules in a magnetic trap at 600 mK 15 . 

The wavefunction for the lowest magnetically trappable Zeeman state of p-CH 2 may be 
written as 

|1 01 ,M S = 1) = a\l 01 ,M s = 1) + Allen, M s = 0) + 7l |l 01 ,M 5 = -1). (33) 

Unlike the ground rotational state of 0-CH2 the spin-rotation and spin-spin interactions 
couple the different and Ms sublevels within the l i rotational manifold of p-CH 2 . Thus, 
the mixing coefficients a\, /3i, and 71 in Eq. (33) are all of similar magnitude, and the cross 
section for spin relaxation in p-CH 2 collisions is given by an expression similar to Eq. (32) 
without the small factor D 2 /AEq 2 (assuming the validity of Born approximation). Thus 
the mechanism of spin relaxation in He + p-CH 2 collisions is similar to that of rotational 
relaxation or rotational depolarization with the corresponding cross sections ranging from 1 
to 100 A 2 (Ref. 31 ). This lack of suppression is responsible for the large magnitude of spin 
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relaxation cross sections for He + p-CB.2 shown in Fig. 8, explaining the large differences 
between the collisional properties of different nuclear spin isomers of CH2 in low-temperature 
collisions with He atoms. 

A simple qualitative explanation for this difference is as follows: Inelastic changes in the 
spin-multiplets correspond to a reorientation of N followed by a recoupling with S, which is 
unaffected by the intermolecular potential which has no magnetic terms. In the case of the 
lowest low-field-seeking Zeeman state of o-CH 2 (|0oo, M$ = +1)), the nonzero N component 
is only a small admixture in the ground-state wavefunction (see Eq. 29) so reorientation 
occurs at a much slower rate. 

D. Isotopic effects 

In the previous section, we have established that the cross section of spin relaxation in 
ground-state 0-CH2 collisions scales as D 2 / AEq 2 with the splitting AEq2 between the ground 
Ooo and the second excited 2o2 rotational levels of 0-CH2. This scaling law is based on a 
simple first Born approximation, which does not take into account the effects of higher-order 
potential couplings and shape resonances (see Fig. 8) on collision dynamics. It is therefore 
desirable to test the accuracy of the scaling prediction, which can be performed by changing 
the rotational constant (and hence the energy gap AE02) of CH2. To this aim, we have 
performed quantum scattering calculations for partially and fully deuterated isotopomers of 
CH2 using the same PES B for all the isotopomers. The rotational and zero-field-splitting 
constants of CH2, CHD, and CD2 used in these calculations are listed in Table II 68 . 

Figure 9 compares the cross sections for spin relaxation in CH 2 , CHD, and CD 2 induced 
by collisions with He at B = 0.1 T. Remarkably, the cross sections exhibit a nearly identical 
dependence on collision energy in the few-partial-wave regime {Eq < 100 mK). This result 
points to the validity of the inverse quadratic scaling of spin relaxation cross sections with 
the rotational splitting (32). Table III lists the ratios of the cross sections for CHD and CD2 
to those of CH 2 calculated at two different collision energies. The ratios are seen to increase 
with a decreasing rotational splitting AE 02 . The deviation from the D 2 / AE^ 2 scaling does 
not exceed 10%, demonstrating the araccuracy of the scaling prediction (32) which is quite 
remarkable given its simplicity. 

In the multiple partial wave regime, the inelastic cross sections in Fig. 9 tend to be dom- 
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inated by shape resonances, which occur due to the centrifugal barriers in the incoming and 
outgoing collision channels 39,70 . These resonances lead to an enhancement of scattering cross 
sections that is correlated with the magnitude of the rotational constant of the metastable 
collision complex. The properties of this collision complex are determined not only by the 
rotational structure of the molecule, but also by the reduced mass and interaction potential 
of the He-molecule complex, leading to pronounced deviations from the simple D 2 /AEq 2 
scaling near scattering resonances. This effect is clearly observed in Fig. 9: In contrast with 
the scaling prediction, the inelastic cross section for CH2 can exceed that of CD2 at selected 
collision energies. We note that similar deviations from the D 2 /AE 2 2 scaling were observed 
experimentally for 3 He + NH (but not 4 He - NH) collisions at 600 K 15 , demonstrating the 
important role of shape resonances in atom-molecule collisions at sub-Kelvin temperatures. 



IV. SUMMARY AND CONCLUSIONS 

We have developed a theoretical methodology for exact quantum mechanical calcula- 
tions of polyatomic molecule collisions in the presence of an external magnetic field. The 
theoretical description is based on a fully uncoupled SF basis set expansion to represent 
the scattering wave function. A system of coupled-channel equations for the radial part 
of the wave function is derived and parametrized by the matrix elements of the molecular 
Hamiltonian and of the atom-molecule interaction PES. Using CH^X) as a general ex- 
ample of a non-linear, asymmetric-top polyatomic molecule, we present the expressions for 
the matrix elements of rotational, centrifugal distortion, spin-rotation, and spin-spin inter- 
actions in terms of the relevant spectroscopic constants 42 (Sec. II and Appendices A and 
B). The matrix elements of the atom- molecule interaction are evaluated by expanding the 
PES in angular basis functions, and we outline the details of this procedure in Sec. II and 
Appendix C. 

We have applied our theoretical framework to explore the dynamics of spin relaxation of 
ground-state CH 2 (A) molecules induced by collisions with He atoms in a magnetic field. The 
scattering calculations were based on two different ab initio PESs developed independently 
of each other: PES A 38 and PES B 31 . We found that, even though the results obtained for 
PES A and B exhibit large quantitative differences, the ratio of elastic to inelastic collision 
rates for He + o-CH 2 collisions greatly exceeds 100 over a wide range of temperatures and 
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magnetic fields regardless of the PESs used to describe the interaction of He with CH2. The 
magnitude of this ratio is thus relatively insensitive to the details of the long-range part of 
the potential, which is described quite differently by the two PESs, as well as to the overall 
"size" of the potential. The large magnitude of the ratio strongly suggests the possibility of 
cryogenic buffer-gas cooling and magnetic trapping of o-CH 2 . 

Our calculations demonstrate that different nuclear spin isomers of CH2 undergo collision- 
induced spin relaxation on widely different timescales. As shown in Fig. 8, the cross section 
for He-induced spin relaxation in p-CH 2 is more than three orders of magnitude larger than 
that for o-CH 2 . These differences arise because the lowest N = rotational state, which 
is more immune to spin relaxation due to the absence of direct intramultiplet spin-rotation 
and spin-spin couplings, is part of the rotational spectrum of o-CH 2 , but not of p-CH 2 . The 
results shown in Fig. 8 thus suggest that it may be possible to isolate the ortho nuclear 
spin isomer of CH 2 by magnetic trapping of the naturally occurring ortho-para mixture. 
Separation of nuclear spin isomers is an active field of research 73 , and our results suggest 
that such separation will occur automatically when a cold ortho-para mixture is confined in 
a magnetic trap in the presence of a cryogenic He gas. 

Using the first Born approximation, we have shown that the cross sections for spin relax- 
ation in the N = rotational state scale as D 2 /AEq 2 with the spin-spin coupling D and the 
rotational splitting AE 02 . This result is fully analogous to the previously derived scaling law 
for diatomic 3 £ molecules 15,72 . We verified that the deviation of the scaling predictions from 
exact CC calculations does not exceed 10% in the s-wave regime (Fig. 9 and Table III). 
The scaling prediction becomes less accurate with increasing collision energy, as the collision 
system enters the multiple partial wave regime (Fig. 9) and shape resonances occur. 

The results of our analysis bear important implications for collisional cooling of poly- 
atomic molecules, opening up the possibility of cryogenic buffer-gas cooling of CH 2 (X) and 
possibly other 38 open-shell polyatomic molecules in a magnetic trap. Experimental real- 
ization of this possibility may open up new directions of research in the study of reaction 
intermediates 71 , large molecule spectroscopy and photochemistry, and may even enable ex- 
ternal field control of complex organic reactions 63 . Moreover, since the magnitudes of the 
elastic-to-inelastic ratios for He + CH 2 plotted in Fig. 8 are similar to those measured 
and calculated for He + NH( 3 £) 15 ' 72 , it is reasonable to assume that collisional properties 
of CH 2 molecules with atoms other than He will be similar to those of the isoelectronic 
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NH( 3 S) molecule 14,15 . If so, it may be possible to use sympathetic cooling with laser-cooled 
alkaline earth or spin-polarized nitrogen atoms 74 ' 75 to create ultracold ensembles of CH2, 
CHD, CD 2 , and possibly other spin-1 and spin-1/2 polyatomic molecules 38 . A more detailed 
theoretical study would be needed to fully explore this possibility, and could be done using 
the formalism presented in this work. 
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APPENDIX A: SPIN-SPIN INTERACTION IN SF FRAME 

Here, we present the derivation of Eq. (6) from the standard MF expression for the 
spin-spin interaction in polyatomic molecules 42 ' 48 ' 49 



(34) 



where the subscripts x, y, and z refer to MF Cartesian components of the spin operators. 
The MF and SF spherical tensor components of S are related by 46 



(35) 



Here and below, the index q (p) is used to denote the MF (SF) components of S 46 . A 
product of two MF spherical tensor operators can thus be expanded as 



Pi,P2 k 



1 1 k 

Qi <h Qi + <h 



1 1 k 

Pi P2 Pi + P2 



Pi+P2,qi+Q2 



(36) 



where we have used a Clebsch-Gordan expansion for the product of two D-functions 46 and 
the quantities in square brackets are Clebsch-Gordan coefficients. Interchanging the order 
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of summation, we obtain 



1 1 k 

qi Q2 qi + Q2 



pi 



i i it 
pi p-pi p 



s (i) s 0) 



Pi "P-PV 



(37) 



where p = p\ + p%. The last sum in this expression is, by definition, a spherical tensor 
product of S with itself 46 , so Eq. (37) takes the form 



c(i)c(i) _ V 



1 1 k 
Qi Q2 qi + <?2 



(38) 



Using the standard relations between the Cartesian and spherical tensor components of 
angular momentum operators 46 , we get three important particular cases of Eq. (38) 66 : 



62 _ 62 
°0 — D z 



E 

fc=0,2 

^=2^=2 x; 

cf Q O q(1) qW 



1 1 k 




l(fc) 
ip ' 



fc=0,2 



fc=0,2 



i i *; 

±1 ±1 ±2 



(fc) 



1 1 k 
±1 =Fl 



(fc) 
p ' 



(39) 



By combining these expressions with the definitions S x = h(S + + S-) and S y = ^(S + — £_) 
and neglecting the terms proportional to S 2 , we obtain Eq. (6) of the main text. 

We note that the following SF tensor expression for the spin-rotation interaction (5) 

2 r 1 

H ST = ^N-S+J2 5(7.-7»)[2^(n)+2^_ a (n)]+^(27,-7.-7 v )2^o(n) [iV<g)S]f, (40) 

p=-2^ - 1 

where 7 = 3(7^ + 7^ + 7z), may be obtained by replacing the second spin operator in Eq. 
(36) with rotational angular momentum iV and retracing the steps outlined above. We 
note, however, that Eq. (5) is simpler and hence more convenient for the practical purpose 
of evaluating the matrix elements in the fully uncoupled SF basis (9). 



APPENDIX B: CENTRIFUGAL DISTORTION MATRIX ELEMENTS 



Since the centrifugal distortion Hamiltonian (4) is spin-independent, its matrix elements 
are diagonal in S and Mg, so we omit these quantum numbers for the rest of this section. In 
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addition, since the MF angular momentum operators do not couple states with different iV 
or Mn (and independent of Mat), we can set N = N' and Mn = M' N in all the expressions 
below. The matrix elements of Eq. (4) in the spherical-top basis \NM^K^) can be evaluated 
by using the standard angular momentum algebra 46 with the result for the first term 



(NM N K N \ - A N N 4 \NM N K' N ) = -A N N 2 (N + 1)H KnK > 



(41) 



the second term 



(NM N K N | - A NK N 2 N 2 | NM N K' N ) 



A 



NK 



N(N+1) 



\±(N, K' N )\ ± (N, K' n ±1)5 Kn , k , n±2 

(42) 



+ 2[N(N+l)-K%]5 KlrK > l 



the third term 
(NM N K N \ - A K N*\NM N K' N ) 
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\ ± (N, K' N )X ± (N, K' N ± 1)\ ± (N, K' N ± 2)A±(iV, K' N ± 3)^, K ' ±4 



2A±(JV, K' N )\±(N, K' N ± l)[2N{N + 1) - (K' ± 2) 2 - K' 2 ]5 Kn ^ n±2 



'2l 



+ {Q[N{N + 1) - K' 2 } 2 - A[N{N + 1) - s, N 
and, finally, the fourth term 

(NM N K N \ - 25 N N 2 (N 2 - N 2 y )\NM N K' N ) = -25 N N(N + 1) 

±A±(J\T, K' N )\±(N, K' N ± l)* WJr±a + {K> 2 - l[N(N + 1) - K*]}5 KlfKi 



(43) 



(44) 



The matrix element of the centrifugal distortion Hamiltonian can be obtained by summing 
Eqs. (41) through (44). 

APPENDIX C: ANGULAR EXPANSION OF THE INTERACTION POTENTIAL 



Here, we derive the expressions for the radial expansion coefficients V^R) based on the 
ab initio PESs V(R, 6, <p) calculated as described in Sec. II. The expansion of the atom- 
molecule interaction PES in spherical harmonics (23) 

Amax A / a \l/2 

VM) + (-inv M (M)], (45) 



A=0 m>0 1 MU \ / 
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may be rewritten in terms of real functions of 9 and <fi as 

Am ax A / o \ 1/2 



V 2A + 1 



<£>\v(&) cos(^), 



(46) 



A=0 /i>0 

where 0x^(6) is a normalized associated Legendre polynomial 46 . A disadvantage of Eq. (46) 
is that the angular basis functions on the right-hand side are not orthonormal. It is more 
convenient to use an alternative expansion in terms of orthonormal angular basis functions 

Amax A |- n "1 1/2 



A=0 m>0 



,1 + <Vo)7T 



cos(/i</>). 



(47) 



The expansion coefficients V\^(R) may be obtained by inverting this expression 



1 


1/2 ^2tt 

/ <^ 


4 / 


.2(1 + ^ )7T. 


Jo 


./o 



sin 6d6V(R, 9, 0)6^(0) cos(^). 



(48) 



The relationship between the coefficients V^i?) and Va^(-R) can be obtained by comparing 
Eqs. (46) and (47). Using Eq. (48), we obtain the final result 

' 2A + 1 \i/2 ,2* 



2tt 



sin 6d6V(R, 9, 0)6 Am (#) cos(//<£). 



(49) 
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TABLE I. Parameters of PES A 



Pairs 


AJeV 


b a /% 1 


C 6 , a /eVa 6 


He-C 


8.67689xl0" 2 


1.03161 


- 1. 30489 xlO 3 


He-H 


24.90525 


1.66630 


1.54137xl0 2 


He-Gi.2 


20.15632 


1.42113 


3.93810xl0 2 


He-G 3 


1.07027xl0 2 


1.82350 


1.99657xl0 2 



TABLE II. Molecular constants of CH2(AT) and its isotopomers used to parametrize the asymptotic 
Hamiltonian (2). All the constants for CH2 are taken from Ref. and expressed in units of cm" 1 . 
The constants for CHD and CD2 are taken from Refs. 43 and 44 . The spin-rotation interaction has a 
negligible effect on spin relaxation cross sections for CHD and CD2, thus the parameters 7^ were 
set to zero. 



Constant 


CH 2 


CHD 


CD 2 


A 


73.06 


55.5243 


37.7869 


D 


8.42 


5.6784 


4.2296 


C 


7.22 


4.9735 


3.6947 


D 


0.778 


0.7567 


0.7765 


E 


0.0399 


0.0461 


0.0406 


lx 


4.46 x 10~ 4 






Iv 


-4.106 x 10~ 3 






lz 


-5.148 x 10" 3 






A N 


3.01 x 10~ 4 


1.97 x 10" 4 


9.295 x 10^ 5 


Ank 


-1.966 x 10" 2 


-7.8 x 10" 3 


-4.964 x 10" 3 


A a - 


1.991049 


1.403 


0.56022 


5n 


1.012 x 10~ 4 


6.19 x 10~ 5 


2.264 x 10~ 5 
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TABLE III. Cross sections for spin relaxation calculated for different isotopomers of CH2 vs the 
rotational energy splitting AE02 (i n units of AE^ICH^)). The cross sections are normalized to the 
He + CH 2 inelastic cross section (<7 ine i = 1.51 x 10~ 4 A 2 for E c = 0.1 cm' 1 and 1.24 x 10 -4 A 2 
for Ec = 0.01 cm -1 ). The values of collision energy (cm -1 ) are given in parentheses. 

Molecule AE 02 cr inel (0.1) cr ine i(0.01) 
CH 2 1 1 1 

CHD 2.19 2.14 2.04 
CD 2 3.98 4.37 3.87 
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Fig. 1. MF coordinate system used to describe the He-CH2 interaction. 




Fig. 2. Contour plots of PES A (left) and PES B (right) for = 0°. 9 and (j) are the spherical polar 
angles of vector R in the MF frame (see Fig. 1 and text for details.) Energies are in cm -1 . 
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Fig. 3. Same as in Fig. 2 for </> = 45°. 



<|> = 90° 

PES A PES B 




30 60 90 120 150 180 30 60 90 120 150 180 
9 / degree 9 / degree 



Fig. 4. Same as in Fig. 2 for <j> = 90°. 



33 




0.5 1 1.5 2 0.5 1 1.5 2 
Magnetic field (T) Magnetic field (T) 



Fig. 5. Lowest Zeeman energy levels of 0-CH2 (left panel) and P-CH2 (right panel). The initial 
states used in scattering calculations are indicated by dashed lines. 
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Fig. 6. Collision energy dependence of elastic and inelastic cross sections for He-CH2 calculated 
using PES A (full lines) and PES B (dashed lines). The magnetic field is 0.1 T. The inset shows 
the thermally averaged ratio of elastic to inelastic collision rates as a function of temperature. 



35 




t — i — r 



-4 -2 -1 

10 10 10 10 10 10 

Collision energy (K) 

Fig. 7. Collision energy dependence of the inelastic cross sections for He-CH2 calculated for several 
values of the magnetic field: 0.01 T (squares), 0.1 T (diamonds), and 1 T (circles). The inset shows 
the temperature dependence of the ratio of elastic to inelastic collision rates for the same values of 
the magnetic field. All the results are obtained using PES B. 
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Fig. 8. Collision energy dependence of the inelastic cross sections for 0CH2 (full lines) and P-CH2 
(dashed lines) in He for B = 0.1 T. The results are obtained using PES B. 
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Fig. 9. Collision energy dependence of the inelastic cross sections for collisions of different methylene 
isotopomers with 3 He: 0-CH2 (diamonds), CD2 (circles) and CHD (squares). The magnetic field is 
0.1 T. The inset shows the temperature dependence of the ratio of elastic to inelastic collision rates 
for the different isotopomers. All the results are obtained using PES B. 
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